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Abstract 


This  report  develops  a  systematic  approach  for  determining  the  acceleration  capability  and 
the  acceleration  properties  of  the  end-effector  of  a  planar  two  degree-of-freedom  manipulator.  The 
acceleration  of  the  end-effector  at  a  given  configuration  of  the  manipulator  is  a  linear  function  of  the 
actuator  torques  and  a  (nonlinear)  quadratic  function  of  the  "joint- velocities".  By  decomposing  the 
functional  relationships  between  the  inputs  (actuator  torques  and  "joint-velocities")  and  the  output 
(acceleration  of  the  end-effector)  into  two  fundamental  mappings,  a  linear  mapping  between  the 
actuator  torque  space  and  the  acceleration  space  of  the  end-effector  and  a  quadratic  (nonlinear) 
mapping  between  the  "joint-velocity"  space  and  the  acceleration  space  of  the  end-effector,  and  by 
deriving  the  properties  of  these  two  mappings,  it  is  possible  to  determine  the  properties  of  all 
acceleration  sets  which  are  the  images  of  the  appropriate  input  sets  under  the  two  fundamental 
mappings.  The  determination  of  the  properties  of  the  quadratic  mapping,  a  key  feature  of  the 
present  work,  allows  us  to  obtain  analytic  expressions  relating  important  acceleration  properties  of 
the  end-effector  to  all  the  manipulator  parameters  and  input  variables  of  interest. 


1  Introduction 


In  this  paper,  we  develop  and  apply  a  systematic  approach  for  studying  the  acceleration  capability  and 
acceleration  properties  of  (a  reference  point  on)  the  end-effector  of  a  planar  two  degree -of-freedom 
manipulator.  The  application  of  the  theory  developed  in  this  paper  to  two  important  problems  which 
arise  in  the  design  of  manipulators  -selection  of  a  manipulator  type  and  determination  of  actuator  sizes  - 
are  described  in  companion  paper  (Desa  and  Kim,  1989).  Acceleration  theory  for  spatial  manipulators  is 
developed  in  a  third  paper  (Kim  and  Desa,  1989). 

An  informal  statement  of  the  acceleration  problem  is  as  follows: 

Consider  the  planar  two  degree-of-freedom  manipulator  shown  schematically  in  Figure  1.  We  are 
interested  in  studying  the  acceleration  of  a  reference  point  P  on  link  2.  ( P  is  typically  a  point  on  the  joint 
axis  of  the  end-effector:  therefore  the  acceleration  of  P  is  often  loosely  referred  to  as  the  acceleration 
of  the  end-effector.)  The  usefulness  of  studying  the  acceleration  of  the  end-effector  of  manipulators  has 
been  discussed  in  (Yoshikawa,  1985),  (Khatib  and  Burdick,  1987)  and  (Graettinger  and  Krogh,  1988)  and 
will  additionally  be  demonstrated  in  (Desa  and  Kim,  1989). 

As  will  be  shown  below,  the  acceleration  capability  of  the  point  P  under  various  conditions  is  best 
described  by  certain  acceleration  sets.  Two  properties  which  arc  used,  in  general,  to  characterize  these 
sets  are  the  maximum  possible  magnitude  of  the  acceleration  of  P  and  the  maximum  magnitude  of  the 
acceleration  of  P  which  is  available  in  all  directions.  The  former  property  is  simply  called  the  maximum 
acceleration  of  P  and  the  latter  the  isotropic  acceleration  of  P  (Khatib  and  Burdick,  1987). 

The  study  of  the  acceleration  properties  of  the  ‘'end-effector”  has  been  a  subject  of  recent  interest 
(Yoshikawa,  1985;  Khatib  and  Burdick,  1987;  Graettinger  and  Krogh,  1988).  It  is  therefore  useful 
to  clearly  state  what  makes  the  problem  of  studying  acceleration  properties  complex  and  how  these 
researchers  have  addressed  this  complexity. 

The  acceleration  of  the  reference  point  P  at  a  given  configuration  (in  the  workspace  of  the  manipulator) 
is  a  linear  function  of  the  actuator  torques  and  a  (nonlinear)  quadratic  function  of  the  rates  of  changes 
of  the  joint-variables  (“joint  velocities").  The  complexity  of  the  “acceleration  problem”  arises  from 
these  quadratic  nonlinearities  in  the  “joint  velocities”.  (Yoshikawa,  1985)  studied  the  acceleration  of  (a 
reference  point  P  on)  the  end-effector  in  connection  with  developing  a  dynamic  manipulability  measure: 
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in  this  study  the  nonlinearities  were  essntially  ignored  since  the  measure  was  estimated  at  zero  “joint 
velocities”.  In  studying  isotropic  acceleration,  (Khatib  and  Burdick,  1987)  dealt  with  the  nonlinearities 
in  a  somewhat  ad-hoc  fashion  by  evaluating  isotropic  acceleration  at  a  “low"  and  a  “high”  joint  velocity 
vector.  (Graettinger  and  Krogh,  1988)  handled  the  nonlinearities  by  posing  the  problem  of  determining 
the  isotropic  acceleration  as  an  optimization  problem. 

In  contrast  to  the  above  approaches,  the  present  paper  demonstrates  how  these  nonlinearities  can 
be  handled  in  an  analytical  manner.  The  fundamental  hypothesis  of  this  paper  is  the  following:  By 
decomposing  the  functional  relationships  between  the  inputs  (actuator  torques  and  joint  variable  rates) 
and  the  output  (acceleration  of  P)  into  two  fundamental  mappings,  a  linear  mapping  between  actuator 
torque  space  and  the  acceleration  space  of  point  P  and  a  quadratic  (nonlinear)  mapping  between  the  “joint 
velocity”  space  and  the  acceleration  space  of  P,  and  by  deriving  the  properties  of  these  two  mappings, 
it  is  possible  to  determine  the  properties  of  all  acceleration  sets  which  are  the  images  of  the  appropriate 
input  sets  under  the  two  fundamental  mappings. 

The  properties  of  linear  mappings  are  well-known.  The  determination  of  the  properties  of  the  quadratic 
mapping  between  the  joint  velocities  and  the  acceleration-space  of  P  is  one  of  the  contributions  of 
this  paper  and  permits  us  to  obtain  exact  analytic  solutions  for  the  isotropic  acceleration  under  various 
conditions. 

In  summary,  the  contributions  of  this  paper  are  the  following: 

1.  Development  of  a  systematic  approach  (stated  in  section  2)  for  defining,  determining  and  charac¬ 
terizing  acceleration  sets. 

2.  Closed-form  analytic  expressions  relating  important  acceleration  properties  of  manipulators  to  all 
the  manipulator  parameters  and  input  variables  (torques,  joint  variable  rates  or  “joint  velocities”) 
of  interest.  (The  only  exception  is  the  maximum  local  acceleration  which  is  specified  in  terms  of 
tight  lower  and  upper  bounds  in  section  6.) 

3.  Necessary  and  sufficient  conditions  for  the  existence  of  isotropic  acceleration.  (Earlier  studies  seem 
to  implicitly  assume  that  isotropic  acceleration  always  exists.)  These  conditions  are  stated  explicitly 
in  terms  of  manipulator  parameters  and  input  variables. 
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4.  Analytical  expressions  for  determining  the  maximum  and  isotropic  acceleration  of  the  end-effector 
at  any  (“local”)  configuration  of  the  manipulator. 

5.  The  theory  treats  nonlinearities  in  an  “exact”  manner  (as  mentioned  above). 

One  consequence  of  2  and  3  above  is  the  development  of  simple  algorithms  (Desa  and  Kim  1989) 
for  sizing  actuators  in  order  to  guarantee  a  specified  isotropic  acceleration.  The  theory  developed  in  this 
paper  is  also  applicable  to  two  degree-of-freedom  manipulators  with  closed-chains  (Desa  and  Kim,  1989). 

The  next  section,  which  describes  our  approach,  also  provides  the  dual  function  of  being  a  “road-map” 
of  the  paper. 


3 


2  Description  of  the  approach 


A  systematic  approach  for  studying  the  acceleration  of  (a  reference  point  P  on)  the  end-effector  based  on 
the  use  of  input-output  mappings  is  as  follows: 

1.  Define  the  input  variables  and  output  variables  of  interest  (subsection  3.1).  The  output  of  interest  is 
the  acceleration  of  the  reference  point  P. 

2.  Define  the  input  sets  of  interest  (subsection  3.1). 

3.  Define  the  input-output  functional  relations.  These  are  obtained  from  the  dynamical  and  kinematical 
equations  of  the  manipulator  (subsection  3.2). 

4.  Define  fundamental  mappings  from  these  functional  relations  (subsection  3.3).  There  are  two 
fundamental  mappings,  a  linear  mapping  and  a  quadratic  mapping. 

5.  Define  the  image  sets  of  the  input  sets  under  the  mappings  obtained  in  set  4  (subsection  3.4).  These 
image  sets  are  the  acceleration  sets  of  interest 

6.  Define  general  properties  which  can  be  used  to  characterize  ("measure”)  acceleration  sets  (subsection 
3.5). 

7.  Determine  the  properties  of  the  mappings  defined  in  step  4  (section  4). 

8.  Determine  the  acceleration  sets  defined  in  step  5  using  the  properties  of  the  mappings  obtained  in 
step  7  (section  4). 

9.  Determine  the  specific  properties  of  the  acceleration  sets  determined  in  step  8  using  the  “measures” 
or  genera]  properties  defined  in  step  6  (section  5). 

10.  Determine  the  local  acceleration  properties  for  any  configuration  q  of  the  manipulator  using  the 
properties  of  the  acceleration  sets  obtained  in  step  9  (section  6). 
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Figure  1:  Schematic  diagram  of  a  planar  two  degree-of-freedom  manipulator 


3  Definition  of  the  acceleration  sets 

3.1  Manipulator  input  and  output  variables 

Consider  a  serial  two  degree-of-freedom  manipulator  with  two  revolute  joints  shown  in  Figure  1.  In  this 
subsection,  we  define  the  link  parameters,  the  input  variables,  the  input  sets,  the  output  variables  and  the 
output  sets  for  a  planar  two  degree-of-freedom  manipulator.  The  manipulator  is  assumed  to  be  rigid  with 
negligible  joint  friction  and  operates  in  a  horizontal  plane  perpendicular  to  the  “gravity  vector”.  (The 
case  of  manipulators  operating  in  gravity  fields  is  relatively  straightforward  and  is  dealt  with  in  (Kim  and 
Desa,  1989).) 

The  link  parameters  necessary  for  describing  the  kinematic  and  dynamic  behavior  of  the  planar  two 
degree-of-freedom  manipulator  (Figure  1)  are  as  follows.  Let  l\  denote  the  length  of  link  1,  a\  the  distance 
from  joint  axis  1  to  the  center  of  mass  of  link  1,  mi  the  mass  of  link  1,  and  /i  the  principal  moment  of 
inertia  of  link  1  with  respect  to  its  center  of  mass  about  an  axis  perpendicular  to  the  plane  of  the  motion. 
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Similarly,  let  k,  <22.  m2,  and  I2  denote  the  corresponding  variables  for  link  2  (see  Figure  1). 

Next,  we  define  the  input  variables,  the  input  constraints  and  the  corresponding  input  sets  of  the  two 
degrce-of-freedom  manipulator.  Let  <71  and  <72  denote  the  generalized  coordinates  of  the  manipulator  (see 
Figure  1),  q\  being  the  joint  variable  at  joint  1  and  <72  the  joint  variable  at  joint  2.  Define 

a  <7t 

Q=  (U) 

<72 

to  be  the  vector  of  joint  variables;  the  corresponding  vector  space  of  all  q  is  called  the  joint  space.  If 


4 iL  <  <7  i  <  <7  iU,  i-  1,2 


denotes  the  constraint  on  joint  variable  i,  then  we  can  define  the  workspace  W  of  a  manipulator  as 

W  =  {q|<7a.  <  <7,  <  qiu,  »  =  1,2}.  (1.3) 

Let  qi  and  72  denote,  respectively,  the  rates  of  change  of  the  joint  variables  <71  and  <72;  ?i  and  qz  will 
be  referred  to  as  joint  variable  rates  for  short.  Define 

.  A  $1 

q=  (1.4) 

<n 

to  be  the  vector  of  the  joint  variable  rates.  If 

I  <7;  l<  iu>,  1=  1,2  (1.5) 

denotes  the  constraints  on  the  joint  variable  rates,  then  we  can  define 


F  =  {q|  I  <7.  I<  <7io,  i=  1,2}  (1.6) 

to  be  the  set  of  all  the  possible  joint  variable  rate  vectors;  graphically  F  can  be  represented  by  (the  interior 
and  boundary  of)  the  rectangle  J\K\hKz  shown  in  Figure  2. 

Let  t\  and  tz  denote  the  actuator  torques,  respectively,  at  joints  1  and  2,  and  define 


A  T\ 
T  = 
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Figure  2:  Set  of  the  joint  variable  rates  of  a  two  degree-of-freedom  manipulator 


to  be  the  actuator  torque  vectors.1  Let 

|  Ti  |<  Tjo,  i-  1,2  O-8) 

denote  the  constraints  on  the  actuator  torques  at  joints  1  and  2.  We  define 

r-{r||n|<rtof  1-1,2}  0-9) 

to  be  the  set  of  the  allowable  actuator  torques;  graphically  T  can  be  represented  by  (the  interior  and 
boundary  of)  the  rectangle  ABCD  in  Figure  3. 

The  vectors  q,  q  and  r  will  be  referred  to  as  the  input  variables  (more  precisely  the  input  variable 
vectors)  of  the  manipulator.  We  will  also  refer  to  the  vector  q  as  a  configuration  of  the  manipulator. 

Let  (xi,  xi)  denote  the  coordinates  of  a  reference  point  P  on  link  2  (see  Figure  1)  in  a  coordinate 
system  fixed  to  the  base  reference  frame  N\  (xi ,  x2)  are  commonly  referred  to  as  task  coordinates.  Define 

(1.10) 

to  be  the  vector  of  task  coordinates;  the  corresponding  vector  space  of  all  xp  is  called  the  task  space. 
‘The  vectors  of  actuator  torques,  joint  variables,  and  joint  variable  rates  denote  column  matrices,  not  physical  vectors. 


Figure  3:  Set  of  the  actuator  torques  of  a  two  degree -of- freedom  manipulator 


The  velocity  xp  and  the  acceleration  xp  of  the  point  P  of  the  manipulator  are,  respectively,  given  by 


xp  = 


x\ 

X-2 


and 


(1.11) 


xp  = 


x\ 

X2 


(1.12) 


The  acceleration  of  P,  xp,  is  the  output  variable  of  interest  in  the  present  work.  The  corresponding  vector 
space  A  of  all  possible  xp  is  called  the  acceleration  space,  expressed  by 


A  =  {x  |  x  G  R2}. 


(1.13) 


3.2  Functional  relations  between  the  inputs  q,  r  and  the  acceleration  xp 

The  next  step  is  to  obtain  the  functional  relations  between  the  acceleration  x*’  and  the  inputs  q  and  r 
for  a  given  configuration  q.  In  this  subsection,  we  show  how  the  necessary  functional  relations  can  be 
obtained  from  the  manipulator  dynamic  equations  and  the  (so-called)  manipulator  Jacobian  relation. 

The  dynamic  behavior  of  the  two  degree-of-freedom  planar  manipulator  in  the  joint  space  can  be 
obtained  using  well-known  methods  (Kane  and  Levinson,  1983;  Kane  and  Levinson,  1985;  Desa  and 
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Roth,  1985)  and  is  described  by  the  following  pair  of  equations: 

d\\'q\  +  dn'qi  -  w  12(^2 +  2q\q2)  =  T\, 

d2\'qi  +  d-nfa  +  W2t<7?  =  T2, 

where  the  coefficients,  dt]  ( ij  =1,2)  and  v,  are  given  in  the  Appendix. 
Defining  the  following  matrix  operators 


D  = 

du 

d\2 

dn 

d22 

W  = 

0 

W12 

W21 

0 

q  = 

q\ 

q2 

r 

1 

{q}2  = 

ii 

<ti  +  2qiq2 

dynamic  equations,(1.14)  and  (1.15),  become 
Dq  + W{q}2  =  r. 


(U4) 

(1.15) 

(1.16) 

(1.17) 

(1.18) 

(1.19) 

(1.20) 


Note  that  equation  (1.20)  is  the  most  general  expression  of  the  dynamics  of  a  two  degree-of-freedom 
planar  manipulator.  The  matrices  D  and  W  standard  for  various  planar  manipulator  types  are  given  in  the 
Appendix.  The  matrix  D  is  the  mass  matrix  of  the  manipulator. 

Since  the  matrix  D  is  always  invertible,  we  can  write  (1.20)  in  a  more  convenient  form  for  our 
purposes  as 

q  =  D-1[r-W{q}2].  (1.21) 

A  crucial  step  in  the  acceleration  analysis  of  a  two  degree-of-freedom  manipulator  is  the  definition 
of  the  matrix  operator  W  and  {q}2,  which  allows  all  the  “non-linearities”  (i.e.  terms  in  the  dynamic 
equations  (1.14)  and  (1.15)  which  are  non-linear  in  the  joint  variable  rates,  q\  and  q 2)  to  be  written  as 
the  product  of  W  and  {q}2.  -  notation  {}2  is  used  to  draw  attention  to  the  fact  that  the  elements  of 
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the  vector  {q}2  are  quadratic  in  the  joint  variable  rates  q\  and  <&.  Note  that  {q}2  is  a  vector  and  should 
not  be  confused  with  the  scalar  q*  which  is  the  square  of  the  magnitude  of  q. 

The  relation  between  the  velocity,  xp,  of  the  point  P,  and  the  joint  variable  rate  vector  q  is  well  known 
(Desa  and  Roth,  1985): 

tf  =  Jq  (1.22) 

where  J  is  a  (2  x  2)  matrix  called  the  manipulator  Jacobian.  The  detailed  expressions  of  the  Jacobian 
matrix  for  various  planar  manipulator  types  are  given  in  the  Appendix. 

To  obtain  the  expression  for  the  acceleration  xp  of  the  point  P,  we  differentiate  equation  (1.22), 


xp  =  jq  +  jq.  (1.23) 

In  the  Appendix,  we  show  that  the  second  terra  in  (1.23),  jq,  can  be  written  in  the  form 

Jq  =  -E{q}2  (1.24) 

where  matrix  E  is  skew-symmetric. 

Substituting  equation  (1.24)  into  (1.23),  we  obtain 

#  =  Jq  -  E{q}2.  (1.25) 

Defining  the  quantities, 

A  =  JD"1,  (1.26) 

B  =  -AW-E,  (1.27) 


it  is  easy  to  verify  that  the  expression  for  the  acceleration  xp  of  the  point  P,  obtained  by  combining 
equation  (1.20)  with  equations  (1.25)  through  (1.27),  is  given  by 

xp  =  Ar  +  B{q}2  (1.28) 

where  A,B  are  configuration  dependent. 

Equation  (1.28)  expresses  the  required  (Input-Output)  functional  relation  between  the  input  variables, 
q  and  r,  and  the  acceleration  xp  of  the  point  P  (the  output  variable)  at  a  given  configuration  q.  It  is 
important  to  note  that  the  definition  of  the  matrix  “operators”  W,  E  and  {q}2  enables  us  to  write  the 
dynamic  equations  in  the  compact  form  (1.28)  which  is  critical  in  the  sequel. 
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33  Mappings 

In  this  subsection,  we  define  two  fundamental  mappings  between  the  input  variables  and  the  acceleration 
\p  of  the  point  P  (the  output  variable). 

It  is  convenient  to  regard  the  functional  relation  (1.28)  as  a  mapping  between  the  input  variables  q 
and  r  and  the  output  variable  xp  for  a  given  configuration  q  of  the  manipulator.  Furthermore,  defining 

=  A  r  (1.29) 

and 


alq 

<*2q 


equation  (1.28)  can  be  written  as 


%p  =  aT  +  aq. 


(1.30) 


0.31) 


The  following  two  simple  and  obvious  relations  are  useful  when  we  define  acceleration'  sets  below: 

xP(q  =  0)  =  qt  =  At  (1.32) 

xp(t  =  0)  =  aq  =  B{q}2.  (1.33) 

It  is  convenient  to  think  of  the  vector  aT  as  the  contribution  of  the  torques  to  the  acceleration  of  the 
reference  point  P,  and  the  vector  aq  as  the  contribution  of  the  joint  variable  rates  to  the  acceleration  of  P. 
The  sum  of  these  two  vectors,  therefore,  gives  us  the  acceleration  of  P  as  expressed  by  equation  (1.31) 
for  a  two  degree-of-freedom  manipulator. 

Equation  (1.29)  can  be  viewed  as  a  linear,  configuration-dependent,  mapping  between  the  torque  vector 
r  and  its  contribution  aT  to  the  acceleration  of  P.  Similarly,  equation  (1.30)  can  be  viewed  as  a  quadratic, 
configuration-dependent,  mapping  between  the  joint  variable  rate  vector  q  and  its  contribution  aq  to  the 
acceleration  of  P.  These  are  the  two  mappings  of  interest  in  this  section. 


3.4  Manipulator  acceleration  sets 

Having  defined  two  fundamental  mappings  of  interest,  we  are  interested  in  the  image  sets  of  the  input 
sets  under  the  mappings  (1.29)  and  (1.30)  at  a  given  configuration  q  of  the  manipulator.  There  are  three 
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image  sets  of  interest. 


3.4.1  Image  set  ST  of  the  actuator  torque  set  T  under  the  linear  mapping 

For  a  given  set  T  of  the  actuator  torques  r  described  by  equation  (1.9),  and  represented  graphically  by  the 
rectangle  ABCD  in  the  r  -  plane  (see  Figure  3),  we  define  the  image  set  ST  of  T  under  the  linear  mapping 
(1.32)  as 

ST  =  {tf’|X'(4  =  0)  =  Ar,rer}.  (1.34) 

(Note  that  ST  lies  in  the  acceleration  plane  A.)  From  equation  (1.32)  and  the  above  definition  (1.34),  we 
see  that  ST  represents  the  set  of  all  possible  accelerations  (the  acceleration  capability  of  the  manipulator) 
when  it  is  at  rest  (q  =  0)  in  any  configuration  q  and  the  actuators  are  turned  on. 

3.4.2  Image  set  Sq  of  the  joint  variable  rate  set  F  under  the  quadratic  mapping 

For  a  given  set  F  of  the  joint  variable  rates  q  described  by  equation  (1.6),  and  represented  graphically 
by  the  rectangle  J\K\JzKi  in  the  q  -  plane  (see  Figure  2),  we  define  the  image  set  Sq  of  F  under  the 
quadratic  mapping  (1.33)  as 

Sq  =  {x'|x'(r  =  0)  =  B{q}2,q  e  F}.  (1.35) 

(Note  that  Sq  lies  in  the  acceleration  plane  A.)  From  equation  (1.33)  and  the  above  definition  (1.35),  we 
see  that  the  image  set  Sq  represents  the  set  of  all  possible  accelerations  (the  acceleration  capability  of  the 
manipulator)  when  the  actuators  are  turned  off  (r  =  0)  in  any  configuration  q. 

3.4.3  State  acceleration  set 

When  a  manipulator  is  in  motion,  the  (dynamic)  state  of  a  manipulator  can  be  specified  by  the  joint 
variables,  (qi,  <72).  and  joint  variable  rates,  (qi,  qi)-  The  state  vector  u  which  characterizes  the  dynamic 
state  of  the  manipulator  is  defined  as  follows: 

(1.36) 


12 


For  a  specified  dynamic  state  of  a  two  degree -of- freedom  manipulator,  the  second  term  of  the  accel¬ 
eration  xp  in  equation  (1.28)  is  a  constant  vector,  which  we  denote  by  k(u)  and  define  as  follows: 

k(u)  = 

Equation  (1.28)  can  then  be  written  as  follows: 

x  =  Ar  +  k.  (1.38) 

For  a  given  dynamic  state  u  of  the  manipulator,  we  define  the  state  acceleration  set,  Su.  as  the  image  set 
of  the  actuator  torque  set  T  under  the  linear  mapping  (1.38): 

5u  =  {xy  =  Ar  +  k,r67}.  (1.39) 

Su  is  therefore  the  set  of  all  possible  accelerations  at  any  given  dynamic  state  u  of  the  manipulator.  Since 
the  dynamic  state  u  of  the  manipulator  essentially  specifies  the  velocity  xp  of  the  point  P  in  (1.11)  in  any 
configuration,  we  can  also  interpret  the  state  acceleration  set  Su  (the  set  of  available  accelerations)  as  the 
acceleration  capability  of  the  manipulator  when  the  manipulator  is  moving  with  the  velocity  xp  in  a  given 
configuration  q. 

3.5  Characterization  of  the  acceleration  sets 

Once  the  acceleration  sets  defined  in  the  previous  section  have  determined,  one  would  like  to  characterize 
them.  In  this  section,  we  define  two  properties  which  are  useful  in  characterizing  acceleration  sets. 

Figure  4  shows  an  acceleration  set  S  in  the  acceleration  plane  x,  and  two  circles  Ci  and  C2.  The  circle 
Cj  of  radius  rj  is  the  smallest  circle  centered  at  the  origin  which  completely  encloses  5.  Its  radius  n 
therefore  represents  the  maximum  (magnitude  of  the)  available  acceleration  in  S.  The  circle  C2  of  radius 
r2  is  the  largest  circle  centered  at  the  origin  which  lies  within  S.  Its  radius  r-i  therefore  represents  the 
largest  (magnitude  of)  acceleration  available  in  all  directions. 

We  define  the  following  two  properties  of  S: 

•  the  maximum  acceleration  of  S:  abi»x(S)  =  r\, 

•  the  isotropic  acceleration  of  S:  Qjjo(S)  =  T2- 


*1 

butf  +  bn[(qi  +<h)2  -  <??] 

*2 

&2i<7i  +  bn[(il\  +  ill)1  -  <7 ?] 

B{q}2.  (1.37) 
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Figure  4:  Characterization  of  an  acceleration  set  in  the  acceleration  plane 


Comments: 

1.  As  will  be  shown,  the  maximum  acceleration  and  isotropic  acceleration  are  two  measures  which 
can  be  readily  extracted  once  the  acceleration  set  is  known. 

2.  The  isotropic  acceleration  (Khatib  and  Burdick,  1987;  Graettinger  and  Krogh,  1988)  is  a  useful 
measure  of  the  acceleration  set,  since  it  is  a  property  which  does  not  depend  on  direction. 

3.  The  average  acceleration  of  the  set  S  cannot  readily  be  extracted  in  closed-form  (or  by  appropriate 
bounds)  from  the  acceleration  set  5.  It  can  however  be  numerically  determined  from  descriptions 
of  the  various  acceleration  sets  given  in  the  next  section.  Also  the  physical  meaning  of  the  average 
acceleration  is  not  clear. 
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4  Determination  of  the  acceleration  sets 

Analytic  expressions  for  the  determination  of  the  three  sets  ST ,  .Sq  and  5U  are  presented,  respectively,  in 
section  4.1,  4.2  and  4.3.  The  determination  of  ST  and  the  state  acceleration  set  Su  follows  directly  from 
well-known  properties  of  linear  mappings  while  the  determination  of  the  set  Sq  requires  the  derivation  of 
the  properties  of  quadratic  mappings  which  are  new 

4.1  Determination  of  the  image  set  5r 

The  set  ST  is  the  image  set  of  the  actuator  torque  set  T  under  the  linear  mapping  (1.32). 

Result  1:  The  image  set  ST  of  the  actuator  torque  set  T  under  the  linear  mapping  (1.32)  is  (the  interior 
and  boundary  of)  the  parallelogram  A'b'C ' D'  in  the  x  -  plane  whose  vertices  A\  b\  C',  and  D'  are 
as  follows: 

A'  :  (aun0  +  a2iV2o,  02\t\0  +  anno), 

B  :  (a\\Ti0  ~  a2iT2o,  a2in0  -  anTio), 

C  :  {-anTio  -  a2\V2o,  -a2\rx0  -  a-nfio), 

D  :  (~a\\Tlo  +  d2\T2o,  -a2\T\o  +  a22T2o),  (1.40) 

where  ciij  (ij=l,2)  are  the  elements  of  the  matrix  A  defined  in  equation  (1.26).  The  centroid  of  the 
parallelogram  A  B' C  D'  is  the  origin  0  of  the  x-plane. 

Result  2:  The  sides  A  B' ,  B'C',  C  D',  and  D  A'  of  the  parallelogram  ST  (Figure  5),  which  comprise  the 
boundary  of  the  set  are  given  by  the  following  equations: 


AB 

:  <Z22*i  -  012X2  =  det(A)ri0, 

(1.41) 

BC' 

:  -a2i*i  +  a\\x2  =  det(A)r2„, 

(1.42) 

CD 

:  a22.i1  -  ai2*2  =  -  det(A)rl0, 

(1.43) 

DA 

:  -a21.il  +  ai  ii2  =  -  det(A)T2*. 

(1.44) 

where  det(A)  is  the  determinant  of  the  (2  x  2)  matrix  A. 
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Figure  5:  Image  set  of  the  linear  mapping  of  a  two  degree -of-freedom  planar  manipulator 


Proof  of  Result  1: 

The  following  are  well  know  properties  of  a  linear  mapping: 

1.  A  line  in  the  r-plane  will  map  into  a  line  in  the  x-plane.  In  particular,  the  line  l\,  with  equation  t\ 
=  0,  maps  into  the  line  4  whose  equation  is 

l\  :  022*1  -  012*2  =  0,  (1 .45) 

and  the  line  lj,  with  equation  T2  =  0,  maps  into  the  line  4  whose  equation  is 

h  ■  -02\*\  +  a\\X2  =  0.  (1.46) 

Both  4  and  4  pass  through  the  origin  (Figure  5). 

2.  Any  line  g\  parallel  to  l\  maps  into  a  line  g[  parallel  to  tv 

3.  Any  line  gi  parallel  to  h  maps  into  a  line  g2  parallel  to  4- 


Regarding  the  rectangle  ABCD  (set  7)  as  a  set  of  lines  parallel  to  li  and  h  one  can  easily  show  the 
well-known  fact  that  the  image  of  ABCD  is  a  parallelogram  A' B' d  D' .  The  vertices  A',  B' ,  d  and  D' 
are  the  images,  respectively,  of  the  vertices  A,  B ,  C  and  D.  Substituting  the  coordinates  of  A(rlo,  7%), 
— n o),  C(-t\0,  -T2o)  and  D(-t\0,  r^)  into  equation  (1.32),  we  obtain  the  coordinates  of  the 
vertices  A' ,  2?',  C‘  and  D'  as  given  in  equation  (1.40).  From  equation  (1.40),  we  see  that  the  vertices 
A'  and  d  are  equidistant  from  the  origin  and  that  the  vertices  B‘  and  D'  are  equidistant  from  the  origin. 
Therefore,  the  origin  of  the  x-plane  is  the  centroid  of  the  parallelogram  a! B' d D' . 

Proof  of  Result  2: 

We  next  need  to  determine  the  equations  of  the  lines  a'b',  B'd,  d D' ,  and  DA,  which  form  the 
boundary  of  the  parallelogram  AB'dD'  in  the  x  -  plane.  A'b'  is  the  image  of  the  line  AB.  whose  equation 
is  r  =  Tia,  to  obtain  the  equation  of  a'b’  ,  substitute  the  equation  of  AB  (r  =  t\0)  into  (1.32)  to  obtain  the 
following  parametric  equations  in  72: 

it  =  anTio  +  ai2T2,  ’  (1.47) 

*2  -  <*2ino  +  aim-  (1.48) 

Eliminating  the  parameter  vi  between  (1.47)  and  (1.48),  we  obtain  the  equation  of  the  line  AB'  in  the  x 
-  plane  as  given  by  equation  (1.41).  In  a  similar  fashion,  we  can  obtain  the  equations  of  the  lines  B'd , 
d D' ,  and  D  A'  as  in  equations  (1.42)  through  (1.44).  Note  that  from  equations  (1.41)  through  (1.44)  we 
see  that  A  B'  is  parallel  to  d D'  and  B'd  is  parallel  to  D  A'  so  that  A  B'dD'  is  indeed  the  parallelogram 
shown  in  Figure  5. 

4.2  Determination  of  the  image  set  Sq 

The  set  Sq  is  the  image  set  of  the  joint  variable  rate  set  F  under  the  mapping  (1.33).  Set  Sq  is  determined 
from  the  following  results 

Results: 

The  set  F  in  the  q  -  plane  is  considered  as  a  family  of  line  segments  passing  through  the  origin.  There 
are  two  such  types  of  line  segments:  those  which  end  on  the  boundaries  J\K\  and  J2K2  parallel  to  the 
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Figure  6:  Image  set  of  the  quadratic  mapping  of  a  two  degree-of-freedom  manipulator 


#2  -  axis,  a  typical  member  of  which  is  the  line  segment  g\  in  Figure  6  (a),  and  those  which  end  on  the 
boundaries  J\Kj  and  J2K\  parallel  tu  .w  q\  -  axis,  a  typical  member  of  which  is  the  line  segment  g2  in 
Figure  6  (a). 


1. 


Every  line  of  the  type  g\  maps  into  a  line  g[  (see  Figure  6  (b))  in  the  x  -  plane,  one  end  of  which 
is  the  origin  and  the  other  end  of  which  lies  on  the  line  segment  f  k  whose  equation  is: 


1  ..  1  ..  (  bn 

biz1  buX2+\  bx2  + 

where  ii  lies  in  the  interval  [Z?u 


=  0 


<fro  +  612WL  +  2<7i<,<72o).  +  biitiH,  ~  2<7i«<72o)]- 


(1.49) 


2.  Every  line  of  the  type  g2  (see  Figure  6  (b))  maps  into  a  line  g2  in  the  x  plane,  one  end  of  which 
is  the  origin  and  the  other  end  of  which  lies  on  the  quadratic  curve  k N1  f  (shown  dashed  between 
k  and  k)  whose  equation  in  the  parametric  (in  ij\)  form  is: 


*i  =  bnifi  +  bndH, +  2^1120), 
h  =  &2iq’?  +  &22(<7i  +  2<?i<72<>). 
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(1.50) 


3.  The  image  set  of  F  is  (the  interior  and  boundary  of)  the  region  Ok  f  k  shown  in  Figure  6  (b),  and 
the  coordinates  of  N1 and  k  are  as  follows: 

k  (  bnfa,  bn&o),  0-51) 

/  (  bn  ~  burfio  +  bnitho  +  fro)2,  (bi\  -bnrf^  +  bniqio  +  qio)2),  (1-52) 

K  (  (bn  -  b n)(}\0  +  bn(q\0  -  q2o)2 ,  (bi\  -  bn)ij\0  +  tm(q\0  -  qio)2)-  (1.53) 


Proof: 

The  quadratic  mapping  is  defined  by  the  following  equation  (1.33): 
x  =  B{q}2 

which  can  be  written  in  the  expanded  form 
*l  =  &n<7i  +  &i2(<72  +  2<?iq2)2, 

X 2  =  62i«i+^2(<?2  +  2qiq2)2-  (1-54) 

The  determination  of  the  set  which  is  the  image  of  the  set  F  (Figure  6  (a))  consists  of  two  steps:2 

1.  Establishment  of  the  properties  of  the  quadratic  mapping,  and 

2.  Determination  of  the  boundary  of  the  image  set  Sq. 

Consider  the  (input)  q  -  plane.  It  is  convenient  to  think  of  this  plane  as  being  generated  by  the 

continuous  family  of  lines  passing  through  the  origin  with  parametric  equation 

' 

q\  =  r 

<  ,  — oo  <  m  <  oo.  (1-55) 

qi~  ml 


Each  value  of  m  gives  us  a  member  of  the  family  of  lines,  a  typical  member  of  which  is  the  line  / 
shown  in  Figure  7.  The  image  (  in  the  x  -  plane  of  the  line  /  is  obtained  by  substituting  (1.55)  into  (1.54) 

J While  the  approach  described  below  is  adequate  for  our  present  purposes,  a  more  basic  approach  to  determining  Sq  is 
described  in  (Kim  and  Desa,  1989). 
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^2 


Figure  7:  Properties  of  the  quadratic  mapping 


and  is  described  by  the  following  parametric  equation, 

3ci  =  [&n  +  b\2(m2  +  2m)]^ , 

x2  =  [ftt  +  &22(m2  +  2m)]r2.  (1.56) 

From  equation  (1.55)  and  (1.56),  one  can  enumerate  the  following  facts: 

Fact  1.  The  image  of  /,  viz.  { ,  is  a  straight  line. 

Fact  2.  The  origin  of  the  q  -  plane  maps  into  the  origin  of  the  x  -  plane. 

Fact  3.  Two  distinct  points  (ft,  ft)  and  (—ft,  -ft)  map  into  the  same  point  of  the  x  -  plane. 

These  results  are  shown  graphically  in  Figure  7. 

Fact  1  follows  from  the  fact  that  (1.56)  is  the  equation  of  a  straight  line  in  the  parameter  t.  Fact  2 
follows  from  the  fact  that  the  point  (0,  0)  in  the  -  plane,  represented  by  the  parameter  t  =  0,  maps 
into  the  point  (0,  0)  in  the  x  -  plane.  If  t  is  the  parameter  corresponding  to  the  point  (ft,  ft)  in  the 
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Figure  8:  Input  6  -  plane 

q-plane,  then  -t  is  the  parameter  of  the  point  (~q\,  -in)  from  (1.55).  From  (1.56),  we  see  that  points 
with  parameters  t  and  -t  will  map  into  the  same  point  in  the  x-plane.  This  proves  Fact  3. 

We  can  therefore  state  the  following  properties  of  the  quadratic  mapping: 

Property  1.  The  image  of  any  line  segment  in  the  q-plane,  one  end  of  which  is  the  origin,  is  a  line 
segment  in  the  x-plane  with  one  end  at  the  origin  of  the  x-plane. 

Property  2.  The  image  of  the  line  /  passing  through  the  origin  in  the  q  -  plane  is  the  half-line  t ,  one 
end  of  which  is  the  origin  (see  Figure  7). 

Property  3.  Given  any  line  l  passing  through  the  origin,  and  the  two  half-planes  p\  and  pi  formed  by  it 
(see  Figure  8),  p\  and  p2  will  have  the  same  image  set  in  the  x  -  plane. 

Property  1  is  a  direct  consequence  of  Facts  1  and  2,  property  2  a  direct  consequence  of  Facts  1,  2  and  3 
and  property  3  follows  from  Fao  3. 

We  now  apply  the  above  properties  to  determine  the  image  set  Sq  of  the  set  F  in  the  q-plane  (See 
Figure  6).  Property  3  tells  us  that  if  we  “bisect”  F  into  two  “half-sets”  with  a  line  passing  through  the 
origin,  then  we  only  need  to  determine  the  image  of  one  of  these  "half-sets”.  The  most  convenient  half-set 


Figure  9:  Input  set  F  in  4  -  plane 

for  our  purposes  is  the  set  K\J\K2  (see  Figure  9),  one  of  the  two  half-sets  formed  by  the  “bisecting”  line 
KXK2. 

In  order  to  use  Property  1,  it  is  convenient  to  view  the  set  K\J\K2  as  a  family  of  line  segments  passing 
through  the  origin;  we  now  need  to  determine  the  image  of  any  line  segment  passing  through  the  origin 
in  this  set.  There  are  two  cases  to  consider:  the  family  of  lines  such  as  gi,  shown  in  Figure  9,  which 
have  one  endpoint  on  the  origin  and  the  other  endpoint  on  the  line  segment  K\J\  and  the  family  of  lines 
such  as  g2.  also  shown  in  Figure  9,  which  have  one  endpoint  on  the  origin  and  the  other  endpoint  on  the 
line  segment  J\K2. 

It  is  convenient  to  decompose  the  half-set  K\J\K2  into  two  subsets  OK\J\  and  OK2J\  as  shown  in 
Figure  10.  Subset  OK\J\  includes  only  the  families  of  lines  such  as  gj  while  subset  OK2J\  includes  only 
the  families  of  lines  such  as  g2.  The  desired  image  set  is  the  union  of  the  images  of  OK\J\  and  OK2J\. 
From  property  1  of  the  quadratic  mapping,  we  know  that  any  line  segment,  such  as  g\  or  g2  of  the  set  F, 
will  map  into  a  line  segment  with  one  endpoint  passing  through  the  origin  of  the  x  -  plane.  To  obtain  the 
other  endpoint  of  the  images  of  the  two  families  of  lines,  we  need  to  find  the  image  of  the  line  segment 
J\K\  of  the  subset  OJ\K\  and  the  image  of  the  line  segment  J\K2  of  the  subset  OJ\K2. 

First,  we  determine  the  image  of  subset  OJ\K\  (10  (b))  by  finding  the  image  of  K\J\.  The  equation 
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of  K\J\  is 


q\=Q\o\  \qi\<qio-  0-57) 

To  obtain  the  image  k  f  of  K\J\,  we  substitute  (1.57)  into  (1.54)  to  obtain 


*i  =  bn<i\o  +  bn(<ii  +  2q\0(n)2, 

*2  =  bixiiXo  +  bz&ifc  +  Tqxoinf.  (158) 

Defining  a  parameter  t  as 

r  =  <?2  +  2q\aq2,  (159) 

equation  (1.58)  can  be  written  in  the  form 

it  =  b\\ij\0  +  bnt,  (1.60) 

*2  =  &2iq?0  +  *22t.  (1-61) 


Eliminating  the  parameter  t  between  the  two  equations,  we  obtain  the  equation  for  the  line  segment 
f  k  as  given  by  equation  (1.49).  (This  proves  the  first  part  of  the  result.)  The  resulting  subset  (Sq)i 
which  is  the  image  of  OK\J\  is  shown  in  Figure  10  (d). 

Next,  we  determine  the  image  of  subset  OJ\Kz  by  finding  the  images  of  KzJ\.  To  obtain  the  image 
fhtk  of  J\Ki,  we  substitute  the  equation  for  J\Kz, 

<72  =  <72o;  Ui  |<  qio  (1-62) 

into  (1.54)  to  obtain  the  parametric  equation  (1.50).  Note  that  (1.50)  represents  the  equation  of  a  quadratic 
curve  in  terms  of  the  parameter  q\.  (This  completes  the  second  part  of  the  result.)  The  resulting  subset 
(Sq)2  which  is  the  image  of  OKiJ\  is  shown  in  Figure  10  (e). 

The  desired  image  set  Sq  of  F  is  the  union  of  (Sq)i  and  (Sq)2  and  is  shown  in  Figure  10  (0- 
Note  that  the  intersection  of  (5q)i  and  (5q)2,  (5q)i  n  (5q)2,  is  not  empty.  Because  Sq  is  the  image 
of  F  under  a  quadratic  mapping,  there  are  points  inside  Sq  which  are  the  image  of  two  distinct  points  in 
K\J\Ki.  In  particular,  any  point  in  the  set  (Sq)i  n  (Sq)2,  will  be  the  image  of  two  points  are  belonging  to 
OK\J\  and  the  other  to  OKiJ\ . 
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Figure  11:  State  acceleration  set  of  a  two  degree-of-freedom  manipulator 

The  images  ,  f  and  k ,  respectively,  of  Wi(0,  qio),  J\(qio,  qio)  and  K\(q\0,  -qio),  are  obtained 
by  substituting  their  (<71 ,  <72)  coordinates  into  equation  (1.54)  to  obtain  the  required  results  (1.51),  (1.52), 
and  (1.53). 


4.3  Determination  of  the  state  acceleration  set  S„ 

The  state  acceleration  Su  corresponding  to  a  state  u  =  (q,q)r  of  the  planar  manipulator  was  defined  by 
equation  (1.39)  and  is  the  image  set  of  the  actuator  torque  set  T  under  the  mapping  (1.38).  We  obtain  the 
following  results  for  the  determination  of  the  state  acceleration  set  Su- 

Result  1:  For  every  element  x(ST)  of  the  image  set  ST,  there  is  a  corresponding  element  x(Su)  of  the 
state  acceleration  set  Su,  given  by 

X(Su)  =  *(ST)  +  k(q,  <*),  (1.63) 
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where 


k(q,  q)  = 


*1 

*ll<7l  +&12(<72  +  2<7l<?2) 

*2 

*>21 9?  +  &22(<?2  + 

B{q}2. 


(1.64) 


Result  2:  The  state  acceleration  set  S„,  corresponding  to  a  state  u  =  (q,  q)r  of  the  planar  two  degree-of- 
freedora  manipulator  is  the  parallelogram  a"b"c" d"  shown  in  Figure  11  obtained  by  translating 
the  set  Sr  by  the  vector  k(q,  q)  in  the  x  -  plane.  The  centroid  of  Su  is  (k\,  h). 


Proof  of  Result  1: 

The  results  1  and  2  are  straightforward. 

From  (1.34),  a  member  x(ST)  of  ST  is  given  by 

x(ST)  =  At.  (1.65) 

From  (1.39),  a  member  x(5u)  of  Su  is  given  by 

x(5u)  =  Ar  +  k  (1.66) 

where  k  is  given  by  equation  (1.64).  Combining  (1.65)  and  (1.66),  we  obtain 

x(Su)  =  x(ST)  +  k,  (1.67) 

which  is  equation  (1.63). 

Proof  of  result  2: 

From  equation  (1.63),  we  see  that  if  we  take  a  vector  x(ST)  of  ST  and  add  the  vector  k  to  it  we  obtain 
the  corresponding  member  x(Su)  of  5U.  So,  if  we  add  the  vector  k  to  every  vector  in  the  set  ST  we  obtain 
the  required  set  5U.  Therefore,  Su  is  the  parallelogram  A"b"C"d”  (Figure  11)  obtained  by  translating  the 
set  ST  (the  parallelogram  A  B  C f  D'  in  Figure  11)  by  the  vector  k.  The  centroid  of  ST  is  x(ST)  =  (0,  0). 
From  (1.67),  we  see  that  the  corresponding  centroid  of  5U  is 

x(Su)  =  0  +  k  =  k.  (1.68) 

This  completes  the  proof  of  Result  2. 
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Figure  12:  Characterization  of  the  image  set  of  the  linear  mapping  for  a  two  de- 
gree-of-freedom  manipulator 


5  Properties  of  the  acceleration  sets 

In  this  section,  we  extract  the  properties  -  defined  in  subsection  3.5  -  of  the  acceleration  sets  ST,  Sq  and 
Su  determined  in  the  previous  section. 

5.1  Properties  of  the  acceleration  set  ST 

We  characterize  the  image  set  ST  of  the  linear  mapping  as  follows: 

Result  1:  The  maximum  acceleration  of  the  acceleration  set  ST  is  denoted  by  an,„(ST)  and  is  given  by 
<Wtfr)  =  max[d«M'),  d{OB))  (1.69) 

where 

d(OA  )  =  \J(a\\T\0  +  axiTio)2  +  feino  +  ariTio)2 

d(OB  )  =  y  {a\\T\o  -  a\2T2o)2  +  (ci2\rlo  -  a22T^)2 
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Result  2:  The  isotropic  acceleration  of  the  acceleration  set  ST  will  be  denoted  by  tfuoO>T)  and  is  given 


ai,o(5T)  =  min[p{A  5  ),  p(SC  )] 


(1.70) 


where 


KAB)  = 


/KB'C')  = 


I  detA  ]  T\0 

V/a12  +  a2i’ 
I  detA  I  T^, 


lan+a2i 


Proof  of  Result  1: 

The  maximum  acceleration  of  Sr  is  the  distance  from  the  origin  to  the  furthest  vertex  of  the  parallel¬ 
ogram  A'b'C ' D'  (see  Figure  12).  Letting  d  ( OA' )  through  d  ( OD ')  denote,  respectively,  the  distances  of 
vertices  A'  through  D'  from  the  origin  O  in  the  x  -  plane,  am«(ST)  is  given  by 


<W(ST)  =  max[  d(OA  ),  d(OB  ),  d(OC  ),  d(OD  ) ). 


(1.71) 


A'  and  C1  are  equidistant  from  the  origin  O.  Also,  B‘  and  D'  are  equidistant  from  the  origin  O.  So, 
<W(ST)  is  given  by 


<W(Sr)  =  max[  d(OA  ),  d(OB  )  ]. 


(1.72) 


Using  (1.32),  the  distance  d(OA')  from  the  origin  O  to  the  point  A' 


d(OA  )  =  (ai\T\0  +  anria)2  +  (a2iT\o  +  azzTTo)2 ■ 


In  exactly  analogous  fashion,  we  obtain 
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The  isotropic  acceleration  of  ST  is  the  shortest  distance  from  the  origin  to  the  sides  of  the  parallelogram 
A'b'C'd'.  Letting  p  ( a'B '),  p  (B'C),  p  (d  D’),  and  p  (D'a')  denote,  respectively,  the  perpendicular 


distances  from  O  to  the  sides  A' B' ,  Bd ,  d D' ,  and  D'a' ,  altQ(ST)  is  given  by 

Oi«,(Sr)  =  mine's'),  p(B'd),  p(C'£>'),  p(D'a')).  (1.75) 

Since  the  origin  O  is  the  centroid  of  the  parallelogram  ST,  parallel  lines  of  the  parallelogram  AB  C'd' 
must  be  equidistant  from  the  origin.  Therefore,  we  can  write  the  following  relations: 

pCAB')  =  p(C'D'),  (1.76) 

p(.B'd )  =  p(D'A).  (1.77) 

Using  (1.76)  and  (1.77),  (1.75)  can  be  written 

Oito(Sr)  =  mm[p(A' B  ),  p(B'd)].  (1.78) 


The  distance  from  a  point  P(Xo,  ya)  to  line  ax  +  by  +  k-Ois  given  by  the  following  well-known  result: 
\axo  +  by0  +  k\ 


dZTP  ‘ 

Using  equation  (1.41),  (1.42)  and  (1.79),  we  obtain 
p(a'b')  = 

Va12  +  ^2 
detA  |  Tfr, 


p(B  C)  = 


+  a21 


(1.79) 


(1.80) 


(1.81) 


Substituting  (1.80)  and  (1.81)  into  equation  (1.78),  we  can  obtain  the  required  result  (1.70)  for  the  isotropic 
acceleration  ftloCS-r). 


5.2  Properties  of  the  acceleration  set  Sq 

We  characterize  the  image  set  Sq  by  the  maximum  acceleration  and  the  maximum  distance  of  any  element 
of  Sq  from  the  two  references  lines  l\  and  h  shown  in  Figure  5. 

Definition  1: 

KtO  =  /(ft,  ft) 

=  \/(b\\<j\  +  b\2<%  +  2bnq\qi)2  +  +  bnife  +  2£22ftft)2  (1-82) 
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^max(x(Sq),  l\ ) 

h 

xi 


Figure  13:  Maximum  distance  from  reference  line  /i  to  a  point  on  the  set  Sq 

Definition  2:  Let  q  denote  the  real  solution  of  the  following  cubic  equation  in  ft: 

[bn$?  +  612WL  +  24i$2o)](£>ii<7i  +  bn<fr) 

+  [£>21  +  ££2($L  +  2tf]fe>)K£>21$l  +  £>22^2)  =  0.  (1.83) 

Definition  3:  Let  p(x(Sq),/ 1)  and  pix(Sq),k)  denote,  respectively,  the  distance  of  any  point  x(Sq)  of  Sq 
from  the  lines  /1  and  fe. 

Pm»*(x(Sq),/i)  =  maxp(x(Sq),/j),  (1.84) 

Pm«(x(Sq),  h)  =  max  p(x(Sq),  h).  (1.85) 

PnuxWSq),  ft),  for  example,  represents  the  distance  of  that  point  of  Sq  furthest  from  l\\  pm%x(x(,S  q),  l\) 
and  Pmw(x(5q),  (2)  are  necessary  for  determining  the  local  isotropic  acceleration  in  subsubsection 
6. 


Definition  4: 


£  a  022611  -  012621  ^ 

022&12  -  Ol2£>22 


(1.86) 
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Definition  5: 


tfifait  $2) 


4 


4 


|  422^1(^1,  fo)  -  fll2*2($l>  fo)  I 

y/M 

I  a-n[b\\tj\  +  fii2(^  +  2q\ fo)]  -  an[b2\ij\  +  bnjtft  +  2<yi <^2)]  I 

V^t  +  °12 

|fl2i*l($l«  fo)  -  QliS2($l,  $2)  I 
>/« 22  + *11 

I  fl2l[fril<7i  +  fci2(<72  +  2<7i<7;)]  -  4n[fc2i<7^  +  £22(92  +  2<?i?2)1  I 

\/*2i+«ii 


(1.87) 


(1.88) 


Result  1:  For  a  general  two  degree -of-freedora  planar  manipulator,  the  maximum  acceleration  of  the 
acceleration  set  Sq  will  be  denoted  by  flmtxCSq)  and  is  given  by 


flmutfq)  =  ntax[/(^i0,  K4\ )  ^2o)i  K$loi  tylo)i  K4loi  ~$2o)]  (1-89) 


where  q\  is  defined  in  (1.83). 

For  the  two  degrec-of-freedom  open-loop  planar  manipulator,  shown  in  Figure  1, 

<2m*x(Sq)  =  Kqiot  Qlo )•  (1-90) 


Result  2:  For  a  general  two  degree -of- freedom  manipulator,  the  maximum  distance  from  an  element  of 
Sq  to  the  reference  lines  /i  and  h  are,  respectively,  given  by 

P««  (  S(SqU)]  (1.91) 

=  max[«7,Wi0,  -$i0),  <r.($",  &<,),  <r,($i<„  $2«),  <r,($l«»  -$2®)L  i  *  1,2.  (1.92) 

where  is  defined  in  (1.86) . 


Proof  of  Result  1: 

The  magnitude  squared  of  the  acceleration  of  a  point  x(Sq)  of  Sq,  denoted  by  J2(5q),  is  given  by 

<?(Sq)  =  ?{q\,  92)  =  Xi(<7i,  ft)  +  *2(91,  to) 

=  (£ll9l  +&1292  + 2^12$!  ft)2 +  (tol9l +^2292  + 2^2291  ft)2-  (1.93) 
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The  maximum  magnitude  squared  of  the  acceleration  for  the  set  Sq,  denoted  by  ajLxCSq).  is  given  by 


aL*(Sq)  =  max/2(<7i,  <fc),  (1.94) 

where  F  is  shown  in  Figure  2  and  specified  by  the  constraints 

\qi  \<  qxo,  (1.95) 

\<n\<<i2o-  (1.96) 

The  maximum  of  Z2,  required  in  equation  (1.94),  will  occur  at  point  q  6  F  which  is  either  inside  F  or 
on  the  boundaries  of  F  where  one  or  both  constraints  might  be  active.  Furthermore,  because  of  property 
3  of  the  quadratic  mapping,  we  only  need  to  look  at  the  boundaries  J\K\  and  J\K2  of  the  half-set  K2J\K\. 
Therefore,  to  obtain  the  maximum  of  (1.82)  under  the  constraints  (1.95)  and  (1.96),  we  should  consider 
the  following  possibilities: 

1.  Neither  of  the  constraints  is  active,  i.e.,  the  max[/2($i,  $2)]  occurs  at  a  point  q  inside  F. 

2.  Constraint  (1.95)  is  active,  i.e.,  max[/2(^i,  <fc)]  occurs  at  a  point  q  lying  on  the  boundary  J\K\  of 
F. 

3.  Constraint  (1.96)  is  active,  i.e.,  maxt/2^,  £2)]  occurs  at  a  point  q  lying  on  the  boundary  J\K2  of 
F. 

4.  Both  constraints  are  active,  i.e.,  max[/2(^i,  qi)]  occurs  at  either  (a)  point  J\  (b)  point  K\  or  (c) 
point  K2.  Since,  by  virtue  of  Fact  3  of  subsubsection  3.1.2,  points  K\  and  K2  have  the  same  image, 
we  only  need  to  consider  either  K\  or  K2\  we  will  choose  K\. 


To  obtain  the  conditions  for  each  one  of  the  above  cases  to  yield  a  maximum,  we  first  differentiate 
/2(<?i  1  Ql)  with  respect  to  q\  and  <72  to  obtain 


+ 


4[b\i<j\  +  &12O72  +  2<7i<72)K&u<7i  +  b\2q2) 

4[&2l9l  +  b 22(^2  +  24?i<72)](^2l4rl  +  &22<?2)> 

(1.97) 

4[hu?i  +  *12(^2  +  2<71<72)]&12(<71  +  $2) 

4[b2l^i  +  ^2(<72  +  lq\<n))b72)(q\  +  <n)- 

(1.98) 
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Now,  wc  consider  each  case. 


Case  1 

To  obtain  the  required  Q  for  the  case  where  both  constraints  are  inactive,  we  set  the  right-hand  side 
of  (1.97)  and  (1.98)  to  zero. 


afi  n  .  di 2  . 

- —  =  0  and  — —  =  0 
dq\  dq2 


(1.99) 


and  obtain 


<7i  =  <72  =  0  (1.100) 

which  actually  corresponds  to  the  minimum  value  of  P{q\ ,  q2),  viz,  zero.  Therefore,  max(/2)  does  not 
occur  at  a  point  q  inside  F. 


Case  2 

Since  constraint  (1.95)  is  active  on  the  boundary  J\K\  of  F,  we  have 

q\  =  q\o  (constant).  (1. 101) 

To  obtain  the  maximum  of  Z2,  we  set  dP/dfc  -  0.  We  therefore  set  the  right-hand  side  of  (1.98)  to 
zero  to  obtain 

q\  +  <72  =  0-  (1.102) 

Combining  (1.101)  and  (1.102),  we  obtain 

<72  =  -q\o  (U03) 

and 

max[/2(<?i ,  q2)]  =  l2(q\0,  -qu)-  (1.104) 


Case  3 

Since  constraint  (1.96)  is  active  on  the  boundary  J\K2  of  F,  we  have 


<72  =  qio  (constant). 


(1.105) 


To  obtain  the  maximum,  wc  set  dpjdqx  =  0.  We  therefore  set  the  right-hand  side  of  (1.98)  to  zero 
and  set  &  =  fo*  to  obtain 

[biitf  +  bi2(.<Ho  +  2(7i<72o)](£ii<?i  +bnq2o) 

+  (&2i<7i  +  +  2<7i0<72)](&21?i  +  bnqio)  -  0.  (1.106) 

Equation  (1.106)  is  a  cubic  in  q\  and  will  therefore  have  three  solutions.  Using  simple  ideas  from  algebraic 
geometry,  we  now  show  that  (1.106)  can  have  at  most  one  real  solution  in  the  region  |c?i|  <  q\0. 

If  l2(qi,  qi)  does  have  a  maximum  4»x  then  the  condition  dP  jdq\  =  0  for  obtaining  4m  is  the 
condition  for  the  quadratic  curve  (1.50)  -the  image  of  J\K2  in  the  x-plane-  and  a  circle  of  radius  4m 
to  have  a  common  tangent  (see  Figure  14).  By  Bezout’s  theorem  (Semple  and  Roth,  1949),  a  quadratic 
curve  and  a  circle  can  have  at  most  two  common  points  of  tangency.  Therefore,  equation  (1.106),  which 
expresses  the  condition  dl2/dq\  =  0,  can  have  at  most  two  real  roots  (one  for  each  point  of  tangency). 
However,  (1.106)  is  a  cubic  in  q\  and  can  therefore  have  either  one  real  root  or  three  real  roots.  Combining 
the  last  two  statements,  we  see  that  (1.106)  can  have  at  most  one  real  root.  (Since  we  are  looking  at 
the  quadratic  curve  in  the  region  |<7i|  <  q\0,  the  real  solution  of  (1.106)  might  lie  outside  the  constraints 
which  simply  means  that  (1.93)  does  not  have  an  extremum  in  the  region  |£i|  <  q\0).  Denoting  the  real 
solution  (1.106)  in  the  region  (<71  (  <  q\a  by  q\  and  using  (1.105),  we  can  write 

ma x[/2(qi,  qi)]  =  f-{q\,  qio )•  (1.107) 

Case  4-a 

When  both  constraints  (1.95)  and  (1.96)  are  active,  and  max[l2(qi ,  <72)]  occurs  at  /i(<7i<„  <720),  then 
max[/2(<7i,  q2)\  =  l2(q\o,  qio)-  (1.108) 

Case  4-b 

When  both  constraints  (1.95)  and  (1.96)  are  active,  and  max[/2(^i,  <72)]  occurs  at  K\(q\0,  -q2o)<  then 

max[/2((7i ,  <72)]  = /2(?i0,  -qio)-  (1.109) 

Therefore,  aman(5q)  is  obtained  as  the  maximum  of  four  quantities  defined  by  equations  (1.104), 
(1.107),  (1.108),  and  (1.109).  This  concludes  the  Proof  of  Result  1. 
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Proof  of  Result  2: 


The  distance  of  any  point  x(Sq)  of  Sq  from  the  line  /„  i=l,  2,  is  given  by 

f*MSq),h)  =  <ri(q\,  <n)  (1.110) 

pW&k)  =  *2(41,  42)  (1-111) 

where  <ti  and  cr2  are  defined  in  equations  (1.87)  and  (1.88).  We  first  wish  to  determine  An«x(x(Sq)./t)  the 
distance  of  l\  from  that  point  of  Sq  furthest  away  from  it  (l\).  pmtx(x(Sq),h)  is  shown  in  Figure  13  for 

given  Sq  and  given  reference  line  l\  and  can  be  defined  as  follows: 


An«(x(Sq),/i)=  max  \o\(q\,  <72)!  (1.112) 

M  (q  en 

where  F  is  shown  in  Figure  2  and  is  specified  by  the  constraints 

I  <71  |<  <7i*,  (1.113) 

I  <72  |<  qio-  (1.114) 

The  maximum  of  <7lt  required  in  equation  (1.112),  will  occur  at  a  point  q  €  F  which  is  either  inside 
F  or  on  the  boundaries  of  F  where  one  or  both  constraints  might  be  active.  Furthermore,  because  of 
property  3  of  the  quadratic  mapping,  we  only  need  to  look  at  the  half-set  K2J1K1  and  its  boundaries  J\K\ 
and  J\Ki.  Therefore,  to  obtain  the  maximum  of  (1.110)  under  the  constraints  (1.113)  and  (1.114),  we 
should  consider  the  following  possibilities. 

1.  Neither  of  the  constraints  is  active,  i.e.,  max<7i(qi,<ft)  occurs  at  a  point  q  inside  F. 

2.  Constraint  (1.113)  is  active,  i.e.,  max  cr\ (q\,  qi)  occurs  at  a  point  q  lying  on  the  boundary  J\K\  of 
F, 

3.  Constraint  (1.114)  is  active,  i.e.,  max o\ (q\,qi)  occurs  at  a  point  q  lying  on  the  boundary  J\Kz  of 
F. 

4.  Both  constraints  are  active,  i.e.,  max  o\ (<71,(72)  occurs  at  either  (a)  point  J\  or  (b)  point  K\  or  (c) 
point  Ki.  Since  K\  and  K2  have  the  same  image  we  only  need  to  consider  K\. 
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To  obtain  the  conditions  for  each  one  of  the  above  cases  to  yield  a  maximum  of  o\ ,  we  first  differentiate 
tfi(<7i?<72)  (equation  (1.87)  with  respect  to  q\  and  $2  to  obtain 

A 

-  2(a22bn  -  anbn)q\  +  2(d22hi2  -  anbn)q\,  (1.115) 

oq\ 

-rr-  =  2(^22612  -  ^12622X^1  +  <?2)-  (1.116) 

dqi 

Now  we  consider  each  case. 


Case  1 

To  obtain  the  required  q,  we  set 

do\  „  .  do\ 

—  =  0  and  —  =  0. 
dq\  dqj 


(1.117) 


We  therefore  set  each  of  the  right-hand  sides  of  (1.115)  and  (1.116)  to  zero  to  obtain 


(71  =$2  =  0  (1.118) 

which  actually  corresponds  to  the  minimum  value  of  tfi(<?i,  $2).  viz,  0.  Therefore,  max(cri($i,  $2))  does 
not  occur  at  a  point  q  inside  F. 


Case  2 

Since  constraint  (1.113)  is  active  on  the  boundary  J\K\  of  F,  we  have 

<7i  =  q\o  (constant).  (1.119) 

To  obtain  the  maximum,  we  set  da\jdq2  =  0.  We  therefore  set  the  right-hand  side  of  (1.116)  to  zero 
to  obtain 

$i+$2  =  0.  (1.120) 

Combining  (1.119)  and  (1.120),  we  obtain 

$2  =  ~q\o  (1121) 

and 

max[<7i($i,  <72)]  =  tfi(<7io,  -<7io)-  (1.122) 
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Case  3 


Since  constraint  (1.114)  is  active  on  the  boundary  J\Ki  of  F,  we  have 

ft  =  <72 «  (constant).  ( 1 . 123) 

To  obtain  the  maximum,  we  set  dcri/dqi  =  0.  We  therefore  set  the  right-hand  side  of  (1.115)  to  zero  and 
set  ft  =  qio  to  obtain  the  linear  equation 


(022*11  -  012*21  )<7l  +  (022*12  ~  012*22)ft  =  0. 

Combining  (1.123)  and  (1.124),  we  obtain  the  solution  of  the  equation  (1.124) 
// 

<?i  =<7  , 


where 

. "  _  022*12  —  012*22  . 
022*11  ~  012*21 

Therefore, 


max[<7 ,(ft,  ft)]  =  <x ,(</",  ft^). 


(1.124) 


(1.125) 


(1.126) 


(1.127) 


Case  4-a 

When  both  constraints  (1.113)  and  (1.114)  are  active,  and  maxoi(ft,  ft)  occurs  at  J\(q\0,  qzo),  then 
max[<7,(ft,  ft)]  =  <Ti(ft<?,  ifio).  (1.128) 


Case  4-b 

When  both  constraints  (1.113)  and  (1.114)  are  active,  and  max<7i(ft,  ft)  occurs  at  K\(qi0,  -qio), 
then 

maxkKft,  ft)]  =  ai(<7w>,  -qio).  (1.129) 

Therefore,  pm»x(x(5q),  l\)  is  obtained  as  the  maximum  of  four  quantities  defined  by  equations  (1.122), 
(1.127),  (1.128),  and  (1.129).  In  exactly  analogous  fashion,  An*x(x(Sq),  h)  is  obtained  as  in  (1.92). 
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Figure  IS:  Characterization  of  the  state  acceleration  set  of  a  two  degrce-of-freedom 
manipulator 

5 3  Properties  of  the  state  acceleration  set  Su 
Definition: 

K  :  centroid  of  the  acceleration  set  in  the  x  -  plane  with  coordinates  Jfci,  £2  given  by  (1.38). 
p(K,  /1)  :  distance  from  point  K  to  the  reference  line  l\. 
p(K ,  I2)  :  distance  from  point  K  to  the  reference  line  h. 

p(AB),  p(A"B"),  :  distance  from  the  origin  to  AB' ,  A"b",  . . .  (see  Figure  16) 

Result  1:  The  maximum  acceleration  corresponding  to  any  dynamic  state  u  of  the  manipulator  is  denoted 
by  am*x(Su)  and  is  given  by 

<W(SU)  =  max[d(04"),  d(OB"),  d(OC"),  dm"))  (U30) 
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where 


d(OA ") 
d(OB) 
d(OC") 
d(QD") 


\J(fl\\T\o  +  a\2T2o  +  ki)2  +  (flixTio  +  022T2O  +  k\)2 
\J(.a\\T\o  -  012^20  + k\)2  +  (a2\T\0  -  <l22T2o  +  k i)2, 

\J \&\\t\o  +  a\2T2o  -  k\P  +  (ai\T\o  +  a22T2o  -  ki)2, 

\J{0\\T\0  -  CL\2T2o  -  kip  +  (a2\Tio  -  022 T2o  -  k\)2. 


Result  2:  The  necessary  and  sufficient  conditions  for  the  existence  of  the  isotropic  acceleration  are  the 
following: 

|  det  A|rl0  -  |*ia22  -  k2an\  >  0,  (1.131) 

|  det  A| Tio  -  |*i<i2i  -  k2a\\\  >  0.  (1.132) 


Result  3:  The  isotropic  acceleration  corresponding  to  any  dynamic  state  u  of  the  manipulator  is  denoted 
by  fliioCSu)  and,  if  conditions  (1.131)  and  (1.132)  are  satisfied,  is  given  by 
|  det(A)|n0  -  | ^22*1  -  a12*2|  |det(A)|T^,  -  |a2i*i  -an*2 1 

V^2+a22  ’  Vau+au 


OvoiSu)  -  min 


(1.133) 


Proof  of  Result  1: 

Let  d(OA ")  through  d(OD")  denote,  respectively,  the  distances  of  vertices  A  through  d"  from  the 
origin  O  in  the  x  -  plane.  Then  a„,„(Su)  is  the  distances  of  the  furthest  vertex  of  the  set  Su  which  is  the 

'V. 

(1.134) 

in  the  x  -  plane  are  given  by 

Xj(A  )  =  xi(A  )  +  k\  =  0)\T]o  +  a\2T2o  +  ki,  (1.135) 

X2  (A")  =  X2(A')  +  k2  =  a2iTi0  +  a22T2<,  +  k2.  (1.136) 


parallelogram  A" B"c” D  " .  Therefore,  amM(Su)  is  given  by 
<W(SU)  =  max[d(OA),  d(OB"),  d(OC"),  d(OD)]. 
Using  (1.40),  the  coordinates  x\(A")  and  'x2(A”)  of  vertex  A" 


The  distance  d(OA  )  from  the  origin  0  to  the  point  A"  is  given  by 


d(OA)  =  ^i\(A")  +  ^{A") 

=  \J (oi  l  Tio  +  a\2T2o  +  k\)2  +  {(i2\T\o  +  a22T2o  +  ki)2. 


(1.137) 
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Figure  16:  Isotropic  acceleration  of  the  state  acceleration  set  of  a  two  degree-of-freedom 
manipulator 


In  exactly  analogous  fashion,  we  obtain 

d(OB  )  =  \J ( auTio  -  anrio  +  k\)2  +  feme  -  a-nTio  +  *i)2,  (1.138) 

d(OC")  =  '/(anrio  +  anrzo  -  ki)2  +  (0MTio  +  a22r2o- ki)2,  (1.139) 

d(OD  )  =  (anrio  -  anrio  -  k\)2  +  ( ai\T\0  -  an -  k\)2.  (1.140) 

Equations  (1.134)  and  (1.137)  through  (1.140)  comprise  result  1. 


Proof  of  Result  2  and  3: 

In  Figure  16,  we  have  shown  two  sets,  ST  and  Su  which  is  obtained  from  ST  by  a  translation  k  = 
(k\ ,  ki)T.  The  centroids  of  ST  and  Sa  arc,  respectively,  by  O  and  K. 
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Using  equations  (1.79),  (1.64),  (1.45),  and  (1.46),  the  distance  from  K  to  the  reference  lines  /i  and  k 
are  given  by 


IKK,  10  - 


(*K,  h)  = 


\j4i  +  an 
1^21^1  -a\\ki\ 


(1.141) 

(1.142) 


\Jd^\+a\\ 

p(K,  l\)  represents  the  perpendicular  distance  between  the  lines  A'b'  and  a"b”  and  also  between  the  lines 
d D'  and  d'D"  (see  Figure  16).  Similarly,  p(AT,  fc)  is  equal  to  the  perpendicular  distance  between  the 
lines  B' d  and  B  'd'  and  also  between  the  lines  DA'  and  D"a"  (see  Figure  16). 

The  state  isotropic  acceleration  aiSO(5u)  is  the  maximum  acceleration  which  is  available  in  all  direc¬ 
tions.  It  is  therefore  equal  to  the  minimum  of  the  distances  from  the  origin  O  (of  the  acceleration  plane) 
to  the  four  sides  of  A"B"d'D"  (the  set  Su). 

Referring  to  Figure  16,  we  can  write  the  following  expression  for  OijoCSu): 


aiso(Su)  =  min[p04"fl"),  p(B"c"),  p(C"d"),  p(D"a")] 


(1.143) 


where  p(A"B")  is  the  (perpendicular)  distance  from  O  to  A"b"  and  similarly  for  p(B"d'),  p(d'D"), 
p(D"a"),  all  assumed  positive  by  definition,  from  the  geometry  of  Figure  16,  we  can  write. 


p(A  B  ),p(C  D  )  =  p(AB)±  p(K,  /i). 


(1.144) 


(Comment:  In  Figure  16,  for  example,  p(A"B ")  =  p(A'B')  +  p(K,li)  and  p(d'D ")  =  p(dD  )  —  p(K,  l\); 
the  correct  choice  of  signs  will  depend  on  the  direction  of  the  translation  but  as  will  be  shown  below  we 
do  not  have  to  worry  about  the  correct  choice  of  signs.) 

Similarly, 


P(B"C"),  p(D"A")  =  p(B'd)  ±  p(K,  l2). 

(The  above  comment  holds  for  (1.145),  too.) 

Combining  equations  (1.143),  (1.144),  and  (1.145),  we  obtain 


(1.145) 


tfiioCSu)  =  m\n[p(A' B  )  ±  p(K,  /,),  p(B'd)±p(K,  l2)l 


(1.146) 
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Since  all  distances  pO  in  the  above  equation  are  positive  by  definition,  we  can  rewrite  the  above  equation 
as 

flboCSu)  =  mine's')  -  p(K,  h),  p(BC)  -  p(K ,  l2)l  (1.147) 

Substituting  equations  (1.80),  (1.81),  (1.141)  and  (1.142)  into  (1.147),  we  obtain  the  required  result 
(1.133). 

Equation  (1.147)  clearly  demonstrates  that  the  isotropic  acceleration  fliW(Su)  for  any  state  u  ^  0  is  less 
than  a„0(ST )  =  min [p(A' B'),p(B'cf)].  In  fact,  if  p(K,  /])  and  p(K,  l2)  are  sufficiently  large  (equivalently, 
the  “nonlinearities”  k\  and  k2  are  sufficiently  “large”),  we  may  not  have  any  isotropic  acceleration.  The 
necessary  and  sufficient  conditions  for  the  existence  of  the  isotropic  acceleration  can  be  obtained  either 
from  (1.147)  or  (1.133).  From  (1.133),  we  obtain  the  following  two  necessary  and  sufficient  conditions 
for  the  existence  of  the  isotropic  acceleration: 

|  det(A)|ri0  -  \k\a-n  -  k2a\2 \  >  0,  (1.148) 

I  det(A)|7^,  -  |*i<z2i  -  Mill  >  0.  (1.149) 

These  are  exactly  the  necessary  and  sufficient  conditions  expressed  in  (1.131)  and  (1.132)  of  result  2. 
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6  Local  acceleration  properties 

At  any  given  (local)  configuration  q  in  the  workspace,  the  following  questions  are  of  theoretical  and 
practical  importance. 

•  What  is  the  magnitude  of  the  maximum  acceleration  at  any  configuration  q  in  the  workspace? 

•  What  is  the  magnitude  of  the  isotropic  acceleration  at  any  configuration  q  in  the  workspace? 

To  answer  both  these  questions,  we  need  to  use  the  properties  of  the  sets  ST,  Sq  Su  developed  in  the 
preceding  subsection. 

Result  1:  The  local  maximum  acceleration  flm«,ioc*i  at  a  given  configuration  q  is  specified  by 

(**max,local)lb  -  «, local  —  (^max, local  )ub  (1.150) 

where  (<w,iocai)ib  is  given  by  (1.130)  with  *i(q,<i)  and  k2(q,(0  evaluated  at  that  joint  variable 
vector  q  which  maximizes  l(q i,  <fi)  in  equation  (1.89). 

(Omax,local)ub  =  <W(Sq)  +  <W(Sr)  (1.151) 

where  Omax(5q)  is  given  by  (1.89)  and  am«x(ST)  is  given  by  (1.69). 

Result  2:  The  local  isotropic  acceleration  local  at  a  given  configuration  q  is  specified  by 

^iio, local 

*  mm[p(A  B  )  -  pmtx(x(Sq), h),  p(BC')  -  pmax(x(5q), /2)1  (1.152) 

where  p(A  B')  and  piB  Cf)  are  given,  respectively,  by  equations  (1.80)  and  (1.81),  and  where 
max(x(5q), /i)  and  max(x(Sq),  h)  are  given  by  equation  (1.92). 

Proof  of  Result  1: 

The  local  maximum  acceleration  Onu,  is  the  maximum  acceleration  over  all  possible  state  acceleration 
sets  Su  at  a  given  position  q  in  the  workspace.  Therefore,  amix  can  be  written  as 

^max.local  =  max(UqgfSu).  (1.153) 
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Figure  17:  Maximum  local  acceleration  of  a  two  degree-of-freedom  manipulator 

It  is  not  possible  to  find  an  exact  analytical  expression  for  a.,,  w.i,  However,  we  can  find  an  upper 
bound  and  lower  bound  which  are  very  good  approximations  to  flmax.iocai. 

Corresponding  to  every  point  P  of  the  set  Sq,  we  have  a  state  acceleration  set  SU(P).  Let  P  be  the 
furthest  point  (from  the  origin)  of  Sq,  and  let  SU(P)  be  the  corresponding  state  acceleration  set,  as  shown 
in  Figure  17.  Also  shown  in  Figure  17  is  the  set  S'U(P')  obtained  by  rotating  the  set  SU(P  )  about  P'  till 
the  longest  diagonal  (A"c"  in  this  case)  of  Su  is  collinear  with  the  line  OP'  joining  the  origin  to  the 
furthest  point  P  of  Sq.  It  is  easily  seen  from  Figure  17,  that  if  vertex  A'  is  the  furthest  vertex  of  SU(P  ) 
from  O ,  then  a  lower  bound  is  given  by 

(<W,ioau)ib  =  d(OA"),  (U54) 


and  an  upper  bound  for  Omax, local  is  given  by 
0*max, local  )ub  =  d(OP  )  +  d(A  P  ), 
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(1.155) 


(*tmax, local  )ub  =  *tmax(5q)  +  flmaaCS-r)-  (1.156) 

In  general,  one  of  the  four  vertices  A"  ,B"  ,C" ,  or  D"  would  be  the  furthest  vertex  of  Su  and  therefore  we 
should  write  (1.154)  as 

(<W,iocj)ib  =  ma x[d(04"),  d(OB"),  d(OC"),  d(OD")].  (1.157) 

Combining  (1.157)  with  equation  (1.137)  through  (1.140),  we  obtain  equation  (1.130).  The  values  of 

and  £2  in  (1.130)  correspond  to  the  furthest  vertex  p  of  S ^  from  the  origin,  i.e.,  to  that  joint  variable 
vector  q  which  maximizes  l(qi,  fa )  in  equation  (1.89).  This  is  simply  a  matter  of  computing  l{q\,  fc) 
at  the  four  vectors  ( q\0 ,  -q\o)T,  (q\,  <ho)T,  (<7io,  qio)T  and  {q\0,  -qzo)r  defined  in  subseaion  3.2.2  and 
determining  which  of  these  four  vectors  maximizes  l(q\,  <72)-  This  completes  the  determination  of  the 
lower  bound  (<w,ioc«i)«»- 

Substituting  for  a,nU(Sq)  and  am*x(ST)  from  equations  (1.89)  and  (1.69),  respectively,  we  obtain 
equation  (1.151).  Thus,  Result  1  is  proved. 

Proof  of  result  2: 

The  local  isotropic  acceleration  is  obtained  in  the  following  steps. 

1.  The  maximum  possible  isotropic  acceleration  is  obtained  when  q  =  0  and  is  equal  to  alu>(Sr )  as 
given  by  equation  (1.70). 

2.  Every  state  acceleration  set  will  have  an  isotropic  acceleration  which  is  less  than  that  given  by 
(1.70)  because  the  “nonlinearities”  effectively  reduce  the  isotropic  acceleration.  The  resulting  state 
isotropic  acceleration  is  a„0(Su)  which  is  given  by  equation  (1.147). 

3.  The  local  isotropic  acceleration  aiJO,ioc*i  is  the  magnitude  of  the  smallest  state  isotropic  acceleration 
at  a  given  local  configuration  q,  i.e. 

^iso, local  =  min  <Zjj0(iSu).  (1.158) 

q  €f 

4.  Using  equation  (1.147)  and  (1.158),  we  can  express  the  local  isotropic  acceleration  <Zj,0,ioc*i  as 

also, local  =  rninmin[p(A  B  )  -  p(K,h),  piBC')  -  p(K,  l2)} 

q  €F 

=  min[min{p(i4  B  )  -  p(K,l\)},  m\n{p(B' C' )  -  p(K\ /2)}].  (1.159) 
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5.  Since  piAB)  and  p{B'd)  are  constants  for  a  given  manipulator  and  given  actuator  constraints, 
(1.159)  can  be  written  as 

tfiio, local  =  min[p(A'B')  -  mdxp(K,l\),  p(B  C')  -  max p(K, I2)].  (1.160) 

where  max(p(K,l\))  is  the  distance  from  the  line  l\  to  the  element  of  Sq  furthest  away  from  /1 
which  we  denoted  in  subsection  3.2.2  by  pm*x(*(Sq),  /j ),  and  max(p(AT,  fe))  is  the  distance  from  the 
line  1%  to  the  element  of  Sq  furthest  away  from  l %  which  we  denoted  by  pmaxCxC-S'q),  h)-  We  can 
therefore  write 

max  p(K,l\)  =  /w(x(Sq),/i)  (1.161) 

maxpiKJi)  =  Pm«(x(Sq),  (2)  (1.162) 

Combining  (1.160),  (1.161)  and  (1.162),  we  obtain  the  required  result  (1.152).  (Note  that  all 
quantities  in  (1.152)  have  been  analytically  determined  earlier.) 


47 


7  Summary  and  conclusions 

In  this  paper,  we  have  developed  a  theory  for  the  acceleration  sets  of  planar  manipulators.  In  particular, 
we  have  accomplished  the  following: 

•  Given  the  kinematical  and  dynamical  equations  of  a  manipulator,  we  have  defined  the  image  set  ST 
corresponding  to  the  set  T  of  actuator  torques,  and  the  image  set  Sq  corresponding  to  the  set  F  of 
the  joint  variable  rates.  We  have  also  defined  the  state  acceleration  set  Su  at  a  specified  point  u  in 
the  state  space. 

•  We  have  determined  the  image  sets,  ST  and  Sq,  and  the  state  acceleration  set  Su. 

•  We  have  characterized  the  image  sets  ST  and  the  state  acceleration  set  S„  by  their  maximum  and 
isotropic  acceleration.  The  image  set  Sq  has  been  also  characterized  by  the  maximum  acceleration. 

•  At  a  configuration  or  position,  q,  in  the  workspace,  we  have  established  two  local  acceleration 
properties:  the  local  maximum  acceleration  and  the  local  isotropic  acceleration.  The  local  maximum 
acceleration  specifies  the  magnitude  of  the  maximum  acceleration  of  (a  reference  point  on)  the 
end-effector,  Hie  local  isotropic  acceleration  specifies  the  magnitude  of  the  maximum  available 
acceleration  of  the  end-effector  in  all  directions. 

We  have,  therefore,  demonstrated  the  hypothesis  which  we  stated  in  the  introduction,  i.e.,  that  the 
analytical  properties  of  acceleration  sets  can  be  determined  from  the  properties  of  the  linear  and  quadratic 
mappings  which  define  them  (the  acceleration  sets).  Furthermore,  the  acceleration  properties  of  interest 
-  especially  the  isotropic  acceleration  -  have  been  determined  in  terms  of  the  manipulator  parameters 
and  the  torque  limits  and  joint  variable  rate  (“joint  velocity”)  limits.  The  stage  has  now  been  set  for 
the  application  of  the  theory  developed  in  this  paper  to  problems  in  the  design  of  manipulators  in  the 
companion  paper  (Desa  and  Kim,  1989). 
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Appendix.  Equations  of  motion  for  the  two  degree-oMreedom  planar  manipulators 

1.  Jacobian  matrix 

The  joint  variable  rate  is  related  to  the  velocity  in  Cartesian  space  by  the  Jacobian  matrix, 

*  =  Jq. 

The  Jacobian  matrix  J  of  the  two  degree -of-freedom  manipulator  shown  in  Figure  1  is  the  following: 

-l\  sin  <71  -  k  sin(<7i  +  <72)  -l2  sin(<7i  +  Ri) 
l\  cos  <71  +  h  cos(<7i  +  <72)  h  cos(<7i  +  <72) 

When  this  relationship  is  differentiated  with  respect  to  the  time,  we  obtain  the  following  equation, 

*  =  Jq  +  )q  =  jq-E{q}2  (1.163) 

where  E  is  the  matrix  which  has  the  following  elements: 

/l  COS  <71  +/2COS(<7i  +  $2)  /2COS(<7i+<72) 

/1  sin  <71  +  l2  sin(<7i  +  q2 )  h  sin(^i  +  q2) 

2.  Dynamic  equation 

The  dynamics  of  the  two-degree-of-frcedom  planar  manipulator  shown  in  Figure  1  is  described  by  the 
following  equation: 

Dq  + V{q}2  =  t  ,  (1.164) 

where  the  components  of  matrices  D  and  V  are  as  follows: 

/ 1  +m\al  + I2  +  m2(al  +  2a2l\  cos<72  +  /i)  h  +  m2(a2  +  a2l\  cos  q2) 
h  +  ^2(^2  +  02k  cos  q2)  I2  +  m2a\ 

0  -m2a2l\  sin  q2 

m2a2l\  sin  q2  0 

and  the  nonlinear  vector  {q}2  is  as  follows: 
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«i}2= 


4? 

(4i  +  42)2  -  4? 


3.  Acceleration  equation 

The  expression  of  the  acceleration  of  the  end-effector  is  as  follows: 
11  =  Ar  +  B{q}2 


where 


A  =  JD"1 
B  =  -AV  -  E 


(1.165) 


(1.166) 

(1.167) 
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